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Abstract. 

We here propose that the large scale superluminal ejections observed in the galactic microquasar GRS 1915-1-105 during 
radio flare events are produced by violent magnetic reconnection episodes in the corona just above the inner edge of the 
magnetized accretion disk that surrounds the central ~ 10 Mq black hole. The process occurs when a large scale magnetic field 
is established by turbulent dynamo in the inner disk region with a ratio between the gas+radiation and the magnetic pressures 
yS 1, implying a magnetic field intensity of ~ 7 x 10* G. During this process, substantial angular momentum is removed 
from the disk by the wind generated by the vertical magnetic flux therefore increasing the disk mass accretion to a value 
near (but below) the critical one (M ~ 10" g s"'). Part of the magnetic energy released by reconnection heats the coronal gas 
(Jc £ 5 X 10* K) that produces a steep, soft X-ray spectrum with luminosity Lx - 10'' erg s"' , in consistency with observations. 
The remaining magnetic energy released goes to accelerate the particles to relativistic velocities (v ~ ~ c, where is the 
Alfven speed) in the reconnection site through first-order Fermi processes. In this context, two possible mechanisms have been 
examined which produce power-law electron distributions N(E) oc with Oe = 5/2, 2, and corresponding synchrotron 

radio power-law spectra with spectral indices which are compatible with that observed during the flares (5,, oc y-o vs.-o.s-j 

Key words, accretion, accretion disks - jets - black hole physics - magnetic fields: reconnection - stars: individual (GRS 
1915+105) 

1. Introduction 

Galactic and extragalactic accretion-powered sources often exhibit quasi-periodic variability and ejection phenomena that may 
offer important clues to the general physical processes in the inner regions of all classes of accreting sources, from young stars 
to X-ray binaries and AGNs. 

The X-ray source GRS 1915-1-105 was the first galactic object to show evidence of superluminal radio ejection (Mirabel & 
Rodriguez 1994, 1999) which is normally interpreted as due to relativistic jets. At a distance of 12.5 kpc and probably with a 
10 solar mass black hole in the center of a binary system, this microquasar offers an excellent laboratory for the investigation 
of black hole accretion and associated jet phenomena revealing, in time scales of minutes to months, a profusion of variable 
emission structures in a large range of wavelengths (Dhawan, Mirabel, & Rodriguez 2000). 

From a detailed analysis of the observed RXTE X-ray variability of this source, Belloni et al. (2000) have distinguished 12 
different modes of variability that they interpreted as transitions between three basic states: a hard spectral-index state (denoted 
by C) with low X-ray luminosity, and two states (denoted by A and B) with softer spectra that were modeled as thermalized 
emission from an accretion disk, with state B being more luminous and hotter than state A. They also found that after a long 
steady hard class, the most frequent class of variability was one characterized by a 30-minute cycle. 

Variability also has been observed in the radio and IR wavebands. Mirabel et al. (1998) reported simultaneous X-ray, IR and 
radio observations of these 30-minute cycles during which the X-ray emission is found to oscillate between states A/B and C. 
The IR/radio emission of the plasma (which is consistent with synchrotron emission from an expanding nuclear jet; see below) 
seems to be associated with the times of recovery from the X-ray dips (see Figs. 2 and 3 from Mirabel et al. 1998), i.e., with the 
times the source is changing from a period in state C to a period in which it rapidly oscillates between states A/B and C. 
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Large scale superluminal radio ejections have been also detected with MERLIN emerging from the source when this was 
at the end of a 3-week extended phase of low/hard X-ray state (also designated the plateau state; Fender et al. 1999; see also 
Fender et al. 2002). Similarly, in a compilation of (VLBA) radio and (RXTE and BATSE) X-ray data, taken during several 
weeks, Dhawan, Mirabel & Rodriguez (2000; hereafter DMR) have distinguished two main states of the system, a plateau and 
a flare state. The plateau state is characterized by a flat radio spectrum {Sy <x v") coming from a compact region of size of a 
few AU, and flux density 1-100 mJy. During this phase, the associated RXTE (2-12 keV) soft X-ray emission is weak, while the 
BATSE (20-100 keV) hard X-rays are strong. On the other hand, during the flare phase, optically thin ejecta are superluminally 
expelled up to thousands of AU with fluxes up to 1 Jy at A\3cm{SyOc y""-*) with rise time < 1 day and which fade after several 
days. The soft X-rays also flare during this phase and exhibit high variability, while the hard X-rays fade for few days before 
recovering (see Figs. 1 and 4 of DMR). The radio imaging over a range of wavelengths (13, 3.6, 2.0, and 0.7 cm) resolves the 
nucleus as a compact jet of length ~ IQXcm AU. This nuclear jet varies in time scales of ~ 30 minutes during minor X-ray/radio 
bursts and reestablishes within ~ 18 hours after a major outburst. So far, GRS 1915 + 105 has been the only system among the 
galactic accreting black holes of stellar mass where both AU-scale steady jets and large scale superluminal ejections have been 
unambiguously observed. 

It is generally believed that the X-ray emission of galactic black holes and active galactic nuclei arises from hot gas accretion 
flow or from an accretion-disk corona. In the first class of models, advection-dominated accretion flow (ADAF) of disk material 
into the black hole (e.g., Narayan, Mahadevan, & Quartaet 1997), as well as, simultaneous advection-dominated inflow-outflow 
(ADIOS) solutions (Blandford & Begelman 1999) have been proposed. In the second class, the presence of magnetic fields 
adds complexity to the solutions but plays an essential role in the disk-coronal coupling, in the energy and angular momentum 
ttansport, and in jet production and acceleration (e.g., Liu, Mineshige & Shibata 2002). 

The AU-scale nuclear jet of GRS 1915+105 developed in the plateau phase has been interpreted as synchrotron emission 
from relativistic electrons (e.g., Falke & Biermann 1999), and its origin has been attributed to disk instabilities in an ADAF-type 
flow by Belloni et al. (1997a, b). According to this model, the 30-minute soft X-ray variabihty observed during the plateau state, 
which involves a transition between states C and A/B, is explained in terms of periodic evacuation and refilling of the inner 
disk region on time scales of seconds as a result of thermal viscosity. This process is believed to occur in the region where the 
radiation pressure dominates the gas pressure and the time required to refill the disk with X-ray emitting hot gas is set by the 
thermal viscous time scale. The observation of radio oscillations with similar period interval supports the idea that the synchrofron 
emission of the nuclear jet is fed by relativistic plasma injected from the inner disk during the dips in the soft X- ray flux (DMR). 

In contrast, the large superluminal radio flares observed in the ~ 500 AU scales cannot be explained by the same viscous disk 
instability model, because they eject an order of magnitude more mass than the AU-scale mini - jet and the recovery time scale 
(of ~ 18 hr) would require according to that model a much larger evacuation radius of disk material well beyond the radiation 
pressure-dominated region where the instability occurs. 

Recently, difi'erent groups have proposed alternative scenarios to explain the variability of the GRS 1915+105 microquasar. 
Livio, Pringle & King (2003; see also King et al. 2004) have suggested that the inner region of the accretion disk would switch 
between two states. In one of them (state A/B), the accretion energy would be dissipated locally to produce the observed disk 
luminosity, and in the other (state C) the accretion energy would be converted into magnetic energy and emitted in the form of 
a relativistic jet (state C). They have attributed the transition between the two states to dynamo generation of a global poloidal 
magnetic field. Similarly, Tagger et al. (2004), also focusing mainly on the 30-minute cycles, have suggested that the alternations 
between the high-soft and the low-hard states and the production of a nuclear relativistic jet could be controlled by the successive 
accumulation and release of poloidal magnetic flux (e.g., by magnetic reconnection) in the innermost parts of the accretion disk. 
They have proposed a specific model to explain the quasi periodic variability in the disk emission based on the development of 
the accretion-ejection instability (Tagger & Pellat 1999). This is a global mode formed by spiral waves that travel back and forth 
in the inner regions of a disk threaded by a poloidal magnetic field with a ratio p between the gas pressure and the magnetic 
pressure close to unity. The spiral waves become unstable by extracting energy and angular momentum from the inner region 
(thus increasing the accretion) and transferring it outward. 

Also, in a contemporaneous work, de Gouveia Dal Pino & Lazarian (2004) have suggested the possibility that the superlu- 
minal large scale radio flares of GRS 1915+105 could be produced by relativistic plasma accelerated during violent magnetic 
reconnection events occurring in the magnetized corona just above the inner regions of the accretion disk. In the present work, 
we try to put together all the pieces of information above in a consistent picture to show that violent magnetic reconnection could 
indeed occur in the innermost parts of the disk and produce the large scale superluminal ejections when p approaches one and 
the mass accretion energy becomes close to the critical value. 

In 12 we draw the potential scenario for flare production and the superlunainal ejecta. Then in %3, we try to quantify it 
by computing the disk and coronal parameters and evaluating the total amount of magnetic energy released during violent 
reconnection in the corona; and in ^4, we evaluate both the electron and the synchrotron power-law spectra that is produced 
by the accelerated particles through first-order Fernai process in the reconnection site. We finaUy conclude (in K5) with a brief 
discussion of the implications of our results. 
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2. A Simple Possible Scenario 

2.1. The BH magnetosphere 

Magnetized accretion disks around rotating (Kerr) black holes (BHs) are frequently invoked to explain the high energy radiation 
and jet production and collimation in quasars and microquasars. The magnetic field lines originally frozen in the disk plasma 
will deposit along with accreting gas onto the BH horizon, therefore developing a magnetosphere around the later (see, e.g., 
Wang, Xiao, & Lei 2002). Near the horizon, no matter how chaotic may be the field threading the disk, the field through the hole 
will become very clean and ordered. If the disk tries to deposit chaotic field on the hole, the field's closed loops will destroy 
themselves on a time scale t ^ Rh/c ~ 10 s Mio, leaving the field ordered. If on the other hand, magnetic pressure from field 
lines threading the disk temporarily push the hole's field into a clumped configuration, it will spring back on this same time scale 
and make itself more uniform (MacDonald et al. 1986). 

The magnetic coupling (MC) of a rotating black hole with its surrounding disk has been recently investigated in some detail 
(see e.g., Blandford 1999, Wang, Lei & Ma 2003, and references therein). The presence of some closed field lines connecting 
the BH with the inner disk edge (in a variant of the Blandford-Znajek process, 1977), can lead to energy and angular momentum 
transfer between the BH and the disk if they are rotating at different angular speeds. Though not a necessary condition, we will 
assume, for simplicity, in the present study that the BH and the inner disk edge are nearly co-rotating so that no significant angular 
momentum and energy transfer is occurring between them. The magnetic structure near the BH horizon will have a force-free 
configuration determined by the MC process and the strength of the field near the hole will be determined by the disk's past 
history. The net flux today must be equal to the total flux that the disk has deposited in its entire past lifetime. These deposits will 
cause the net flux through the hole to fluctuate stochastically. The rms field will grow with time until it becomes so strong that 
it pushes its way back into the accretion disk and slows down accretion. The resulting maximum strength is likely to be such as 
to produce equilibrium between gas ram pressure and magnetic pressure in the innermost parts of the disk (Wang, Lei and Ma 

2003) , though Parker-Rayleigh-Taylor instabilities may reduce this strength (MacDonald et al. 1986). 

2.2. The accretion disk and corona around the BH 

In the standard model for an accretion disk around a black hole, gravitational energy is continually converted into kinetic, thermal, 
and magnetic energy by the joint action of viscosity, and magnetic field Une stretching and reconnection. The viscosity can be 
mainly turbulent or magnetic, or both (Shakura-Sunyaev 1973). 

In the process of angular momentum transport in the disks, shear is expected to convert an existing poloidal magnetic field 
into a toroidal component. Then, probably some of the toroidal flux will be used to generate new poloidal fields by raising 
magnetic loops in a dynamo process. Most of the energy will be either dissipated within the disk that will then emit thermahzed 
radiation, or in the low density corona above and below the disk that will radiate thermally or by Compton scattering, or else the 
energy will be released in the form of a bulk flow, as a jet. In the last case, it is possible that essentially all the accretion energy 
of the inner disk is lost in the jet (see, Blandford & Payne 1982, Livio et al. 2003 and references therein). 

Let us assume as in previous works (e.g., Merloni 2003, Livio et al. 2003, Tagger et al. 2004, de Gouveia Dal Pino & Lazarian 

2004) that the inner parts of the accretion disk alternate between two states which are controlled by changes in the global magnetic 
field. Most of the magnetic field generated in a dynamo process in which the turbulence is driven by, e.g., the magneto-rotational 
instabiUty (Balbus & Hawley 1991), is a tangled small-scale magnetic field within the disk. As remarked by Livio et al. 2003, 
as long as the magnetic field remains as such most of the dissipation will occur locally within the disk that will be radiative and 
very luminous. However, if a large scale poloidal magnetic field is generated then it is possible that an energetic outflow can be 
driven. 

Numerical simulations of the dynamo process in accretion disks (e.g., Kudoh et al. 2002, Shibata 2003), are still unable 
to compute the generation of a global poloidal magnetic field, so that for the time being we must rely on physical estimations 
rather than on detailed modeling (Livio, Pringle and King 2003). Tagger et al. (2004), for example, have suggested that a gradual 
increase in the magnetic flux in the inner regions of the accretion disk could be achieved by the combined action of the turbulent 
dynamo associated, e.g., with the magneto-rotational instability (which would create a vertical flux of opposite polarities in the 
irmer and outer disk regions), followed by advection of the magnetic flux with the gas to the inner disk regions. Though magnetic 
diflrisivity could interrupt this magnetic flux accumulation, there is at least some observational evidence for it in the case of the 
magnetic field structure of our own galactic disk that shows a strong vertical field in the inner tens of pes of the galactic center 
which is almost three orders of magnitude more intense than the horizontal magnetic fields measured elsewhere in the disk. 
Livio (1997; see also Livio et al. 2003, King et al. 2004), in addition, assuming a turbulent dynamo process with inverse cascade 
and reconnection of magnetic loops above and below the disk, have demonstrated that (provided that the distribution nil) of the 
lengths of the magnetic loops satisfies a relation «(/) oc Z^^/^) large scale poloidal fields can be generated such that Lj ~ Lace, 
where Lj is the jet power and Lace is the disk accretion luminosity. 

Considering the aspects above, we assume here that during the plateau state that precedes a radio flare, a large scale poloidal 
field is progressively built in the disk by a turbulent dynamo process. The action of buoyancy forces will also make the disk 
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Fig. 1. Schematic drawing of the magnetic field geometry in the inner disk region at Rx- The acceleration occurs in the magnetic 
reconnection site at the Y type neutral zone (see the text). 

unstable against the Parker-Rayleigh-Taylor instability and horizontal magnetic field lines will raise from the disk forming large 
scale loops in the rarified hot corona. We further assume that once the dynamo process establishes a global poloidal field over 
a substantial region of the disk, it will be able to maintain that field for a period of time. The length of time the source spends 
in this state is hard to estimate without a complete physical model, however observations indicate that the system remains in the 
plateau state for prolonged periods of time. 

The vertical field flux will give rise to a wind that will remove angular momentum from the disk, therefore significantly 
increasing the accretion rate (possibly at a rate greater than the rate due to the disk viscosity). Also, with the accumulation 
(trap) of vertical flux in the inner regions the ratio between the gas+radiation pressure to the magnetic field pressure (J3) will 
soon decrease to one. Under these conditions, events of reconnection of magnetic field lines with opposite polarization will be 
inevitable, and in the innermost regions this process may become eventually very violent when enough magnetic energy is stored 
in the corona. We show below that this occurs when p ~ I and the accretion rate approaches the critical value. We then argue 
that this could explain the large scale radio flares in GRS 1915+105 with the ejection of the observed relativistic blobs from the 
innermost regions. 

2.3. Magnetic reconnection site in the inner disk-corona region 

Considering the assumptions above, the resulting quaUtative structure of the magnetosphere of the hole and the accretion disk 
must be like that shown in Figure 1 (which was adapted from Macdonald & Thorne 1982, and Macdonald et al. 1986). In the inner 
disk region (of radius Rx), there is a site which may be appropriate for violent magnetic reconnection. Surfaces of null poloidal 
field lines mediate the geometry of the open field lines anchored into the BH horizon with the opened lines of the disk wind 
and those connecting the disk with the BH horizon. Labeled as 'T netral zone" in Figure 1, these magnetic null surfaces begin 
or end on Y points. Across each null surface, the poloidal field suffers a sharp reversal of direction. According to the Ampere's 
law, large electric currents must flow out of the plane shown in Figure 1, along the null surfaces, and in the presence of finite 
electric resistivity, dissipation of these currents will lead to reconnection of the oppositely directed field Unes (e.g., Biskamp 
1997, Lazarian & Vishniac 1999). We here investigate the possibility that the magnetic energy released by reconnection near the 
Y point region (see also Figure 2) is able to accelerate the plasma to relativistic velocities and produce the synchrotron radio 
flares (see below). 

The release of magnetic energy by reconnection is believed to be the most promising mechanism for production of solar 
flares. Both direct and indirect astrophysical evidence show that the nature seems to know the way of producing fast 
reconnection in spite of the theoretical difficulties involved (see, e.g., Lazarian, Vishniac & Cho 2004). Y-type neutral 
zones are also present, for example, in eruptive magnetic arcades of the solar corona, like those associated to two-ribbon solar 
flares (see, e.g., Takasaki et al. 2004, Sakajiri et al. 2004). In these cases, a possible interpretation of the observations is that a 
current sheet forms between an emerging magnetic flux and the pre-existing magnetic flux of opposite polarization and magnetic 
reconnection in this zone will cause the formation of a filament/plasmoid. This will temporally inhibit further reconnection and 
establish a temporally stable magnetic configuration. During this stage, there will be considerable magnetic energy storage due 
to the continuing converging motion of the two fluxes. Eventually an eruption develops (for example, from sheared or twisted 
arcades lying underneath the magnetic null point), the filament/plasmoid is then ejected and a large quantity of magnetic energy 
is released (see, e.g., Takasaki et al. 2004). 
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Lazarian & Vishniac (1999) proposed a mechanism at which flares can naturally emerge as a result of an instabil- 
ity mediated by turbulence. In their model the reconnection rate increases due to magnetic field wandering caused by 
turbulence. Fluxes of non-parallel magnetic field lines reconnect very slowly if the turbulent velocities are much smaller 
than the Alfven velocity. However, outflows caused by reconnection increase the level of turbulence and therefore the 
reconnection rate. As a result of such a positive feedback a flare develops. 

Whatever is the particular mechanism, the energy released by magnetic reconnection is expected to ultimately drive violent 
outward motions in the surrounding plasma that can produce copious amounts of solar cosmic rays via direct acceleration along 
the open magnetic field lines or via a first order Fermi process (e.g., Reames 1995). Similar mechanisms have been recently 
proposed to accelerate the plasma in magnetic reconnection zones of young born accretion-induced collapse pulsars to produce 
ultra-high energy cosmic rays (de Gouveia Dal Pino & Lazarian 2000, 2001). 

The several stages described above for magnetic reconnection and solar flare production may be also operating in our ac- 
cretion disk-corona-BH system, but since in the present study it is not crucial to follow the detailed evolution of the magnetic 
configuration and of the features that develop during the reconnection process, we will consider a simpler scenario. According 
to the discussion in the previous sections, we argue that whenever becomes small enough and the accretion rate attains values 
near the critical one in the inner disk region, then the lines of opposite polarization near the Y-zone will be pressed together by 
the accreted material sufficiently rapidly as to undergo violent reconnection events. Energy stored in the magnetic field will be 
released suddenly to heat the plasma and accelerate charged particles. Some of the high speed particles will travel downward to 
the BH or to the disk to the foot points of the magnetic-arches, yielding bright X-ray emission, others will spew outward giving 
rise to the relativistic radio ejecta (see below). ' 

We notice that the merging of the closed field lines of the loops and the open field lines in the disk further out from the 
inner radius can also lead to reconnection, however, we find that the amount of magnetic energy released in these regions is not 
sufficient to accelerate the particles to produce the observed radio luminosities in the large flares (see below). 

For reconnection to occur in the inner disk region Rx, near the Y point, the poloidal flux from the disk and the flux anchored 
in the central magnetic structure surrounding the BH must have opposite signs, like in Figure 1 . As pointed out by Tagger et 
al. (2004), this suggests that large ejections with violent flares will occur only when the system is in that configuration. Violent 
reconnection will not take place in the periods when both fluxes have the same sign. This dichotomy is observed for example in 
the interaction between the terrestrial and the solar wind magnetic field. Now, the magnetic flux trapped in the BH horizon is a 
sum of all magnetic flux accreted during the whole history of the source. If the magnetic flux in the disk is either produced by a 
turbulent dynamo in the disk itself or advected from a companion star, its sign is expected to change from time to time, although 
the characteristic time scale of this variation is difficult to compute in the absence of a detailed modelling of the dynamo in the 
disk or in the companion star On the long term, a sufficient number of reconnection events might in turn cancel the flux trapped 
in the BH horizon and change its sign as well. In the case of the solar dynamo, it is not clear yet why the field reversals take 
place every 22 year. In the case of GRO 1915-1-105, Tagger et al. (2004) have suggested that this long term evolution of the field 
configuration could be of the order of one to few years so that part of the time the field fluxes in the inner disk edge and in the 
hole would be antiparallel and the rest of the time parallel. 



' In the scenario just described, the high accretion rate may also cause the compression and shear of magnetic loops that may exist in the 
lower parts of the inner disk corona, under the Y-zone dome. These could also undergo violent reconnection and even produce plasmoids that 
could eventually erupt and provoke particle acceleration as well, like in the Sun (see discussion in section 4). 
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As remarked above, observations of the Sun indicate that the amplification, dissipation and reconnection of the magnetic field 
lines is likely to occur not steadily, but sporadically, involving violent flares followed by more gradual recoveries to the pre-flare 
magnetic field configuration with no obvious periodicity. Similar behavior, i.e., violent flares followed by more slowly recovering 
of the pre-flare state is also detected in the observed superluminal radio ejecta of GRS 1915+105, thus suggesting that magnetic 
reconnection could be also associated to the production of relativistic blob ejections that rapidly move to the ~ 100 AU scales. 
In the next section, we try to quantify the ideas just described. 

We notice that after reconnection, the partial destruction of the magnetic flux in the inner disk will make /? to increase and 
make the disk resumes the less magnetized condition of the plateau state with most of the energy being dissipated locally within 
the disk instead of in the outflow. 



3. Building tlie IVIodel for tlie Superluminal Ejections During the Flares 

We adopt the field geometry of the magnetized accretion disk/corona as described in the previous section (Fig. 1). 

We could first try to make an order-of-magnitude estimate of the amount of magnetic energy that could be extracted from the 
inner disk region by simply assuming that it is approximately equal to the disk accretion energy rate GM*M//?x- For a BH with 
mass Mi, - lOM© = Mio and Schwartzschild radius /J* = 2GMi,/c^ = 2.96 x 10'' Mio cm, this results 

Wb - 10^^ ergs-' Mio Mig/?^;, (1) 

where we have assumed that the inner radius of the accretion disk (Rx) corresponds approximately to the last stable orbit of the 
BH, Rx ^ 3Ri, ^ 10^ cm, and that the accretion rate M < Msdd, where MEdd = 1 .9 x 10''* Mio g s ' is the Eddington critical 
accretion rate (Shakura & Sunyaev 1973), M19 is M in units of 10^^ g s-\ and Rx,i is Rx in units of lO' cm. This estimate of Wb 
can be compared to the observed X-ray luminosity of GRS 1915+105 (~ 10^^ erg s"'). 

Now, in order to evaluate more precisely the amount of magnetic energy that can be extracted through violent reconnection 
from the Y neutral zone to produce superluminal ejecta, we need first to evaluate the physical conditions both in the magnetized 
disk and corona. 



3.1. Disk and Coronal Parameters 

To evaluate the disk quantities at Rx, we may employ the standard disk model (Shakura & Sunyaev 1973; Frank et al. 1992). 
In the inner disk region (named the X-region), the radiation pressure is larger than the thermal pressure and the disk will be 
radiation pressure-dominated. In this case, 

Td ^ 1.2x10^ K a~\^'^ mI(^^ R'^J,"^ (2) 
and 

Ud ^ 2.8 X 10'^ cm-3 1 M'l'^ M'^ r][1j q-^\^ (3) 

where, q - [I - (/?*//?x)'^^]'^'* = 0.82, and are the disk gas number density and temperature, respectively, ao.i is the 
Shakura- Sunyaev viscous coefficient in units of 0.1 (e.g., Frank et al. 1992), and M19 is the disk accretion rate in units of lO'' g 



The width H of the disk in the radiation pressure-dominated region is given by 

^ ^0.7M,9R-x^jqt,,2 (4) 
Rx 

As remarked in the previous section, the violent magnetic reconnection events that may originate large scale relativistic 
ejections are expected to occur when the large scale magnetic field in the inner disk region, which is generated by dynamo 
process, attains an intensity as large as B^/Sn ^ Pd, where Pd is the disk pressure which is dominated by the radiation pressure, 
i.e., whenjS ^ 1, giving 

Bd ^ 6.9 X 10^ G a~^^ mJ^'* R'^.j'^ (5) 

Recently, Liu, Mineshige, & Shibata (2002), have proposed a simple model to quantify the parameters of a magnetized disk- 
corona system. Assuming, like in the solar corona, that gas evaporation at the foot point of a magnetic flux tube quickly builds up 
the density of the corona to a certain value, and the tube radiates away the heating due to magnetic reconnection through Compton 
scattering, they derived the following coronal quantities (as a function of disk quantities) for a radiation pressure-dominated disk 

- 7.6 X 10'^ cm-3 a-\'"'j3;"' M-f'^^ M'^ R^f' if" (6) 
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and 

^ 4.9 X 10« K alf'^p]"^ M'f M^l 4'/;^ I'J' (7) 

where ric and Tc are the coronal gas number density and temperature, respectively, ySi is the value of yS in the disk in units of 1, and 
is the scale height of the Y neutral zone in the corona in units of 10^ cm. Although the quantities above have been computed 
for the purposes of the present model only in the inner disk radius Rx, they are actually valid for any radius of the disk-corona 
within the radiation pressure-dominated region. We note that, according to the Liu et al.'s model, smaller values of yS would result 
more energy dissipation through the corona and larger coronal temperature and density, as indicated by Eqs. (6) and (7). 



3.2. Rate of Magnetic Energy Release in the Y-type neutral zone 

Let us assume that the poloidal magnetic field that raises in the corona just above the inner disk region (Bx) is of the order of 
the local disk magnetic field, Bj (eq. 5), that has been produced and pushed to the inner disk region as described in the previous 
section (see also, e.g., Arons 1993, de Gouveia Dal Pino & Lazarian 2001). ^ 

The rate of magnetic energy that can be extracted from the Y-zone in the corona (above and below the disk) through recon- 
nection is Wb - {B\l%n) ^va i^^Rx) Lx, where ^ = Vrec/vA is the reconnection efficiency factor, Vrec is the reconnection velocity, 
va = Z?x/(47rnemp)'^^ is the coronal Alfven speed, nip is the hydrogen mass, and Lx is the length of the reconnection region 
which according to Figure 3 is Lx/ARx ~ va/^va ~ ^"^ For the conditions here investigated we find that va - c. Substitution 
of these relations into the equation above gives 

B^ 

Wb ^ -^ci4nRx)ARx (8) 
8;r 

In order to estimate ARx (or Lx), we need to model the reconnection process itself. Classical reconnection schemes use to 
provide too slow reconnection rates. However, there is ample observational evidence that reconnection in Astrophysics may take 
place at very large rates which are comparable with the Alfven speed (see e.g., Dere 1996). In the framework of solar flares, 
considerable progress has been made in recent years towards the determination of more realistic reconnection rates during flare 
events based mostly on the separation speed of Ha flare ribbons (e.g., Qiu et al. 2004; Jing et al. 2005; Sakajiri et al. 2004, 
Takasaki et al. 2004). These observations indicate that the reconnection speed can be as high as few tenths of the Alfven velocity 
(Takasaki et al. 2004). Since in the conditions we deal with in the present analysis this speed approaches c, the corresponding 
reconnection rate can be rather large. 

As mentioned earlier, in the turbulent reconnection model proposed by Lazarian & Vishniac (1999, henceforth LV99; 2000, 
see also Lazarian 2005 for an extensive review on this subject) wandering of magnetic field lines was suggested as the ultimate 
cause of fast reconnection. An alternative model (Bisckamp, Schwarz & Drake, 1997, Shay et al. 2003) attempts to explain 
fast reconnection appealing to plasma properties, such as anomalous resistivity (Parker 1979). These models do not ex- 
clude each other. For instance, while dealing with anomalous resistivity in most astrophysical situations one has also to 
deal with the effects of turbulence. At the same time the reconnection rates in the LV99 model may get an initial boost 
from the anomalous resistive effects. This boost may be important when the initial turbulent velocities are much smaller 
than the Alfven speed. In our present analysis, we focus our attention on regions of strong magnetic fields, like those 
around the BH in the X-region during nearly critical accretion rate. For those regions the resistivity should be anomalous 
and the fast reconnection should be guaranteed (with Vrec ~ va). 

Following LV99, we can estimate the width of the current sheet for which the resistivity should be anomalous: 

cAB 

ARx = (9) 

47mcZev,h,c 

where AB ^ 2Bx denotes the change of the magnetic field across the reconnection region, and v,h,c = (5/MZi)i/2 the thermal 
velocity of the ions of charge Ze in the corona 

v,u,c = 2.65 X 10«cm s"' ^ /3f"' M^^"^' M^l'^ Rlf' (10) 
Substitution of Eqs. (5), (6), and (10) into Eq. (9) results 

I Mix 



- 0.0034 Z-i a^f Vf^" m'^^I'^' R-f"'' 4'/'^ (1 1) 



R^ 



^ We notice that this value is a little larger than the one estimated by assuming, e.g., that the ejected gas takes the bulk of the accretion power, 
i.e., Lj ~ Lace (where Lj ^ BIr\v^I2, Lace - Mv^, and is the Keplerian toroidal velocity in the disk; see also Livio et al. 2003). 
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Fig. 3. Schematic representation of a reconnection region. 

Now, if we plug the result above and Eq. (5) into Eq. (8), we get 

. 2.4xlo3Vgs-ao1«? (^) (12) 
and the corresponding reconnection time is 

tree - ^ - 3.3 X 10-4 S-1 Rx,7 (13) 

which indicates that the release of magnetic energy is very fast. 

The analysis of Eqs. (11) and (12) shows that the adoption of a somewhat smaller viscosity parameter a (which would imply 
a more quiescent disk), by an order of magnitude, for instance, would require a disk accretion rate only few times larger in order 
to keep the values of ARx/Rx and Wb unaltered. 

Eq. (12) gives the total expected amount of magnetic energy rate released by fast reconnection in the Y-zone during the flare 
state of GRS 1 9 1 5 + 1 05 . Part of this energy will heat the gas that will produce X-ray luminosity, and the remaining will accelerate 
the charged particles to relativistic velocities producing a violent radio ejecta. But what is the fraction that goes to each part, and 
what is the potential acceleration mechanism in the reconnection site that can produce the superluminal radio ejections? 

The evaluation of the exact fraction of the magnetic energy that goes to each of the processes above is a task out of the 
scope of the present work, as it depends on the detailed modelling of the conditions of the disk-corona and of the efiiciency of 
each process. However, assuming that most of the magnetic energy released through recoimection in the corona goes to heat the 
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plasma, Liu et al. (2002) have computed the average value of the Compton parameter (y) for an optically thin corona with an 
original distribution of thermal electrons, y = ^^HcCTtI, which was weighted over all the disk radii. According to their results, 
for a disk near the critical accretion (M ~ MEdd), y < 0-8, which implies a soft (and steep) X-ray spectrum, just as required by 
GRS 1915+105 observations during the flare state, with a spectral index cKr-.t = (| + f) ^ i > This result is qualitatively 
consistent with the present model, as it predicts, for a magnetized disk-corona a predominantly soft X-ray spectrum during nearly 
critical accretion. 

Since observations do indicate that both, heating and acceleration of the plasma to relativistic velocities are occurring, as 
suggested by the observed disk X-ray luminosity of the order of ~ 1 0^^ erg s" ' , and by the radio flares which in turn suggest a jet 
mechanical power also of the order of ~ 10^^ erg s"^ at least in the nuclear region (see, e.g, Falcke & Biermann 1999), we will 
here assume that a substantial fraction of Wb (Eq. 12) goes to accelerate the particles. 

Let us now try to answer the second part of the question above, i.e., what could be the acceleration mechanism in the 
reconnection region able to produce the relativistic ejecta and the radio flare? 



4. The Acceleration Mechanism during the Flare State 

As stressed by de Gouveia Dal Pino & Lazarian (2001), the particular mechanism of particle acceleration during reconnection 
events is still unclear in spite of numerous attempts to solve the problem (see LaRosa et al 1996, Litvinenko 1996). Cosmic rays 
from the Sun confirm that the process is sufficiently efficient in spite of the apparent theoretical difficulties for its explanation. 

Particles could, in principle, be accelerated at once to the relativistic energies by the large induced electric field within the 
reconnection region (e.g., Bruhwiller & Zweibel 1992, Haswell, Tajima, & Sakai 1992, Litvinenko 1996; de Gouveia Dal Pino 
& Lazarian 2000, 2001), but one-shot acceleration mechanism can make it diflftculty to explain how particles will afterwards 
produce a power-law particle and synchrotron spectrum. Another possibility is that particles are accelerated by a Fermi process 
in the reconnection site. In this regard, we discuss below two viable sorts of this process. 

A schematic representation of a reconnection region is shown in Figure 3. The upper and lower parts of the magnetic flux 
move towards each other with the velocity v, . As a result, charged particles rays in the upper part of the reconnection zone "see" 
the lower part of the magnetic flux to approach them with a velocity 2vrec- In this case an acceleration process analogous to the 
first-order Fermi acceleration of cosmic rays in magnetized shocks may take place (e.g., Longair 1992) and we can write 

AE 2Vrec 

— = —^cose (14) 
E c 

For a particle probabihty distribution p{6) = 2 sin cos 6d6, the average energy gain per crossing of the reconnection region is ^ 

AE Vrec ^Vrec 

( — ) = -^ I 2cos^esinede= -— (15) 
E c Jq 3 c 

where 6 is the pitch angle between the particle velocity and the magnetic field. The round trip of the particles between the upper 
and lower magnetic fluxes will produce an increment of energy twice as large. Thus the ratio of the particle energy, E, after such 
a round trip to its original energy, Eq, will be 

5=A = i 8^ (16) 
Eq 3 c 

Now, assuming that the escape probability of a particle from the acceleration zone is similar to the one computed for a shock 
front, 4vrec/c, then probabihty to stay within the reconnection region is simply 

4Vrer 

P=l —, (17) 

c 

and the electron particle spectrum is obtained from (see e.g., Longair 1992) 

NiE)dE ~ E-^^'^dE oc E-^l^dE (18) 

This equation indicates that the particle spectrum produced by first-order Fermi acceleration within the reconnection site is steeper 
than that normally produced in shocks (for which a similar simplified treatment that ignores losses and non-linearities provides 
A'(£') oc E^^), so that the two processes may be, in principle, distinguished. Now, the Synchrotron radiation from electrons with 
a power-law distribution N{E) oc E~p, has a power-law frequency radio spectrum oc v~"\ where the spectral index is given 



^ We note that since the particle Larmour radius around the magnetic field in the reconnection zone, vg = ynieC^/iZeBx) ^ 9 x 10"* cm 
73.5 Bxj « ARx (where 73.5 = y/3.5 is the Lorentz factor for electrons with velocities ~ 0.92 c), the particles will remain confined time 
enough in the reconnection zone to be efficiently accelerated through the Fermi mechanism before escaping from the system. 
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byar - ip - l)/2 - 0.75. Though crude, this estimated value is consistent with the observed radio spectrum during the flares 
(ar ^ 0.6). 

Finally, the particle acceleration mechanism discussed above does not disregard the possibility that the relativistic fluid may 
be also produced behind shocks. As in the Sun, plasmoids formed by reconnection of the field lines lines near the Y-zone may 
violently erupt and cause the formation of a shock front. This raises the alternative possibility of a first-order Fermi acceleration 
process of the ionized gas behind the shock front and the development of a power-law spectrum with particle spectral index of 
the order of two (NiE) oc E^^) as remarked above and, as such, a corresponding synchrotron spectrum with a spectral index also 
of the order of the observed one, Sy oc v~^-^. 

After a radio flare, while the system is slowly recovering its original magnetic field configuration, the plateau phase must 
resume. During this recovering phase, the magnetic field intensity in the X-region will be smaller and also the mass accretion 
rate. In this state, the X-ray spectrum is observed to be hard and flat. As we have seen in 13.2, the model of Lin et al. (2002) 
predicts a harder X-ray spectrum when the accretion rate in the magnetized disk is very sub-Eddington for which a spectral index 
ar-x < 1.2 should be expected. 

5. Conclusions and Discussion 

We have here investigated the origin of the large scale superluminal ejections observed in the galactic microquasar GRS 
1915+105 during radio flare events and proposed a very simple model in which they could be due to violent magnetic reconnec- 
tion episodes in the corona just above the inner edge of the magnetized accretion disk that surrounds the central ~ IOMq black 
hole. In this region, when a large scale magnetic field is built by a turbulent dynamo process, and there is approximate balance 
between the disk magnetic and gas+radiation pressures (J3 - 1), which implies a strong magnetic field intensity - 7 x 10^ G, 
reconnection between the magnetic lines that arise from the disk with those of the magnetosphere of the BH may become violent. 
During this regime, the vertical magnetic flux generates a wind that will substantially remove angular momentum from the disk 
therefore increasing the disk mass accretion to a value near the critical one (withM ~ 10''' gs-'). 

Whenever the system reaches the condition above and the magnetic field lines of the BH magnetosphere and those of the 
disk have opposite polarization, it may undergo an episode of violent reconnection that lasts for very short time ~ 3 x lO"'' s. A 
substantial part of the total magnetic energy released by reconnection (~ 2 x 10^^ erg s"') heats the coronal gas to ^ 5 X 10^ 
K that produces a steep, soft X-ray spectrum with the expected luminosity from observations. Assuming that the remaining 
released magnetic energy goes to accelerate the particles to relativistic velocities (v ~ va ~ c, where va is the Alfven speed) in 
the reconnection site through first-order Fermi processes, we have found that a power-law electron distribution N(E) oc E^^^^-^^, 
and a synchrotron radio power- law spectrum (5y oc y-o.75,-0.5-) produced with radio spectral indices which are comparable to 
the one observed during the flares. 

After each violent reconnection event, the system must resume a sub-critical disk mass accretion with smaller magnetic field 
intensity, and a slower reconnection rate, thus recovering the more quiescent plateau state, as described in fl to 4. 

Although until the present, GRS 1915 + 105 has been the only system among the galactic accreting black holes where large 
scale (> 100 AU) superluminal ejections have been observed (maybe because this system is the only one that meets the necessary 
conditions as investigated here), the magnetized-disk-corona model here discussed could be applicable, in principle, also to other 
X-ray binary sources, or even to the microquasar extragalactic counterparts, i.e., quasars with evidence of superluminal ejection. 
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